Directed flow at midrapidity in heavy-ion collisions 
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It was recently shown that fluctuations in the initial geometry of a heavy ion collision generally 
result in a dipole asymmetry of the distribution of outgoing particles. This asymmetry, unlike the 
usual directed flow, is expected to be present at a wide range of rapidity — including midrapidity. 
First evidence of this phenomenon can be seen in recent two-particle correlation data by STAR, 
providing the last element necessary to quantitatively describe long-range dihadron correlations. 
We extract differential directed flow from these data and propose a new direct measurement. 

PACS numbers: 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

Analyses of correlations between particles emitted in 
heavy-ion collisions reveal azimuthal structure that is not 
seen in proton-proton or deuteron-gold collisions. These 
unique correlations are present even if the particles are 
separated by a large interval in rapidity. The largest 
component, known as elliptic flow, is one of the early ob- 
servations at RHIC [1]. More recently, it was noticed 
that more detailed features of the correlation pattern 
known as "ridge" and "shoulder" can also be explained 
by a new phenomenon, dubbed "triangular flow" , which 
comes from a hydrodynamic response to fluctuations in 
the initial geometry [2] ■ An analysis Q of recent experi- 
mental data [4J shows that all the correlations observed 
at large relative pseudorapidities are likely to originate 
from collective flow (and global momentum conserva- 
tion). Contributions from event- by-event fluctuations, 
which in the past have been neglected, are key to this 
understanding. 

Teaney and Yan [f| have shown that such fluctuations 
in the initial geometry are expected to create a new type 
of directed flow in addition to elliptic and triangular flow. 
Fluctuations break the symmetry of the initial density 
profile, and as a result there is, in general, one direction 
where the profile is steepest. This effect can be quantified 
as a dipole asymmetry in the initial density [5(: 
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where the averages in the right-hand side are taken over 
the initial transverse entropy density profile, and (r, 4>) is 
a polar coordinate system around the center of the distri- 
bution, chosen such that {re 1 ^) — 0. If one chooses E\ to 
be positive, then $i generally corresponds to the steepest 
direction for a smooth profile, and s± is the magnitude 
of the dipole asymmetry. In general, e± will differ from 
— even at midrapidity — due to fluctuations. 

This dipole asymmetry, followed by hydrodynamic ex- 
pansion, creates a specific type of flow. Recall that par- 
ticles with large transverse momentum pt are created 
where the fluid velocity is largest, and their momentum 



is parallel to the fluid velocity |fj. The fluid velocity 
scales like the gradients of the initial density profile, and 
so high-p t particles are likely to be emitted along the 
steepest gradient, i.e., with azimuth $i. Just as the 
quadrupole asymmetry in the distribution of outgoing 
particles is quantified by the well-studied elliptic flow co- 
efficient V2, the second Fourier harmonic with respect to 
the event plane, this asymmetry can be characterized by 
the directed flow measured with respect to $1 
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where <p p is the azimuthal angle of the outgoing parti- 
cle momentum and the brackets indicate an average over 
these particles in a given collision event. 

While vi is positive for high-j>t particles, the condition 
that the total net transverse momentum of the fluid ap- 
proximately vanishes implies in turn that v\ is negative 
for \ow-pt particles, which results in a specific pattern for 
the pt dependence of v\ @. Although $1 cannot be mea- 
sured experimentally, the correlation of every particle to 
<&! induces a correlation between pairs of particles, which 
can be measured even if $1 is not known. 

Unlike the usual directed flow, which is odd in rapid- 
ity 0-HJ , this new type of directed flow is expected to de- 
pend little on rapidity, since it is created by fluctuations 
in the initial geometry. Hints of this phenomenon can al- 
ready be seen in existing data, and it enables a complete 
quantitative understanding of long-range dihadron corre- 
lations (in conjunction with the previously studied ellip- 
tic flow, quadrangular flow 0, triangular flow, and mo- 
mentum conservation) |3| . In this article, we extract dif- 
ferential directed flow from recent experimental data [4j , 
and we propose a method to measure it directly. 



II. DIHADRON CORRELATIONS 

The STAR collaboration has released an analysis of di- 
hadron correlations at different angles with respect to the 
event plane [4] . They have measured the distribution of 
the relative azimuthal angle A</> between a high-p t trigger 
particle and an associated particle in various pt bins. 
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Along with the contribution from flow, this measure- 
ment also contains a correlation induced by the global 
conservation of transverse momentum. The overall prob- 
ability of detecting a pair of particles with transverse 

momenta p[ and p[ a ^ can be written in terms of the 
one-particle distribution dN/cPpt, under fairly general 
assumptions, as (Toj : 
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where the sum is over all particles in the event and the 
brackets average over all events in the centrality bin un- 
der consideration, and p t = \/pt- The term in parenthe- 
ses represents the preference for back-to-back emission 
due to the near-vanishing total transverse momentum. 

The one-particle distribution is typically characterized 
in terms of Fourier components 
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Plugging this into Eq. ([3]), the first Fourier harmonic of 
this two-particle correlation, (cos A/fi) , can then be writ- 
ten, to first order, as a sum of contributions from mo- 
mentum conservation and from directed flow Q 
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Here, and vf 1 are the directed flows of the associated 
and trigger particles, and higher order terms are negligi- 
ble when each of these is much smaller than 1 . Explicitly, 
then, the momentum conservation correlation is 
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We assume that these are the only contributions at large 
relative pseudorapidity A77. In practice, we use STAR 
data projected onto |A^| > 0.7. 

In order to isolate the flow contribution, one must 
quantitatively estimate the effect of momentum conser- 
vation. We estimate the denominator of Eq. ([6]) in the 
following way: we use the total pion and kaon multi- 
plicities measured by BRAHMS in central collisions [ll[ , 
and we rescale them to the 20%-60% centrality range 
assuming that the correlation ^ scales with centrality 
approximately like 1/-/V part , which increases the correla- 
tion by a factor ~ 3.6. We calculate (pf ) for pions, kaons 
and nucleons assuming exponential mt spectra, and using 
the inverse slopes of m t spectra measured by PHENIX at 
mid-rapidity in the centrality range 40%-50% We 
further estimate that the effective inverse slopes, aver- 
aged over all rapidities, are 10% smaller than the inverse 
slopes at mid-rapidity, based on the observed decrease of 
(pt) versus rapidity We thus obtain 
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FIG. 1. (Color online) Differential vi of charged hadrons ex- 
tracted from STAR correlation data 0]. Circles (squares) cor- 
respond to trigger particles with p^*' = 3-4 (4-6) GeV. Closed 
(open) symbols correspond to assumed values of =0.1 
(0.13). Error bars are statistical only. The value of pt on the 
horizontal axis is the average value in the corresponding bin. 
For p t > 0.5 GeV/c, this value is obtained using p t spectra 
measured by PHENIX [TJ| in the centrality bin 30-40%. For 
the lowest bin, it is obtained by assuming an exponential ra t 
spectrum. For sake of clarity, (pt) has been shifted by —0.05 
(0.05) for open (closed) symbols. 



Subtracting the contribution of momentum conserva- 
tion from the measured cos A(f>, we extract v^v^ using 
Eq. (JSJ) - Since v\ is defined as a correlation with a di- 
rection $1 which is itself uncorrelated with the reaction 
plane [5|, one expects that this quantity is independent 
of the orientation of the trigger particle. This is indeed 
what is observed by analyzing STAR results Q. In or- 
der to increase the statistics, we average the measured 
(cosA0) over all orientations. We thus obtain 10 val- 
ues of (cos A4>) , corresponding to 5 intervals in p t below 
3 GeV/c for the associated particle, and 2 intervals in p t 
for the trigger particle, namely, = 3-4, 4-6 GeV. 



Assuming some value of , we can extract v\ a> from 

the product v^v^. There is some arbitrariness in this 
procedure, but v\ is expected to be a smooth function 
of pt, so that it should have comparable values in the 
intervals 2-3 and 3-4 GeV. Any value of in the range 
0.10-0.13 gives reasonable results for both sets of trigger 
particles, as shown in Fig. [TJ In particular, we obtain 

similar curves for v[ a ^ from both sets of trigger particles, 
which supports the validity of Eq. ^ . 

Teaney and Yan Q used ideal hydrodynamics to make 
predictions for this V\, and, although the values given 
were not intended as precise, realistic predictions, the 
generic properties seen in Fig. [T] match these expecta- 
tions. The variation of v\ with p t is as predicted below 
2 GeV (see Fig. 12 of Ref. [|). In particular, the ab- 
sence of net transverse momentum implies that the sum 
of ptViijPt) over all particles vanishes, which explains why 
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v\ changes sign near 0.8 GeV/c. Above 2 GeV, v\ devi- 
ates from the linear rise predicted by ideal hydrodynam- 
ics and saturates or decreases, as expected, in the same 
way as elliptic flow [l4| , due to viscous corrections and / or 
the onset of hard physics not captured by hydrodynam- 
ics. The magnitude of Vi is also roughly as expected: 
the initial dipole asymmetry E\ defined by Eq. ([1]) is pre- 
dicted to be approximately 4 times smaller than the ini- 
tial eccentricity £2 for Au-Au collisions [5J. In the mo- 
mentum range 1.5 < p t < 2 GeV, hydrodynamics pre- 
dicts Vi/si ~ V2/s2- Since V2 — 0.2 in this momentum 
range for the same centrality window (l4| . one expects 
vi ~ 0.05, in good agreement with the value in Fig. [TJ 

The estimates in Fig.[T]have sizable errors. First, there 
are unnecessary statistical errors (shown in the figure), 
because they are extracted from correlations involving 
a trigger particle with high transverse momentum, and 
there are few such particles. Uncertainty in (p t ) for each 
bin is also present. Most important, however, is the large 
systematic error coming from the estimation of the mo- 
mentum conservation correlation. We estimate that our 
value for the denominator in Eq. ([5]) is only accurate to 
~ ±20%. This can change the extracted value of v\ by 
as much as 50%. 

Fortunately, these errors can be significantly reduced 
by carrying out a dedicated analysis. 



respect to midrapidity. By construction, this results in 
a measured v\ that, on average, vanishes at midrapidity 
and is antisymmetric in y. 

In this paper, we study directed flow created by fluc- 
tuations in the initial geometry. These fluctuations are 
expected to depend weakly on rapidity [l8|]. Thus we 
choose a weight that depends on p t only. The choice 

w(pt)= Pt -^, (9) 

where angular brackets denote an average over parti- 
cles in the detector acceptance, satisfies the condition 
(wpt) = 0, so that momentum conservation does not bias 
the analysis. It also corresponds to the expected behav- 
ior of v\{pt) in ideal hydrodynamics in the limit of low 
freeze-out temperature for massless particles 

Instead of the event-plane method, one can use any of 
its variants such as the scalar-product method All 
methods should give the same result provided that one 
only correlates particles, or subevents, separated by a gap 
in pseudorapidity in order to remove nonflow effects. 



IV. CONCLUSIONS 



III. NEW METHOD OF ANALYSIS 

Directed flow is most often analyzed using an event- 
plane method, analogous to the one often used for elliptic 
flow [15l.ll6l]. In every event, one defines the directed flow 
"event plan^B' ipEP,i by 

Q COS IpEP.l = J]] Wj COS 

Q sin ipEP.i = J]] Wj sin cf>j (8) 

where the sum is over particles detected in the event, <pj 
are the azimuthal angles of outgoing particles, and Wj is 
a weight, and Q > 0. Our new method is based on a 
new choice for the weight Wj . The weight of a particle is 
generally a function of its transverse momentum p t and 
rapidity y. The desired conditions are that 

• Correlations from momentum conservation should 
not bias the flow analysis. Since these correla- 
tions are proportional to pt , this gives the condition 
(w(pt,y)pt) — 0, where angular brackets denote an 
average over the detector acceptance. 

• Particles with larger v\ should be given more 
weight, the optimal choice being w(j>t,y) oc 

For smooth initial conditions, the symmetry between tar- 
get and projectile implies that v\ is an odd function of y. 
A standard choice is w(pt, y) = y, which satisfies both re- 
quirements if the detector acceptance is symmetric with 



We have shown that the first Fourier harmonic, 
(cos A0) , of the two-particle correlation measured by 
STAR shows the first evidence for a sizable directed flow 
originating from fluctuations of the initial geometry. This 
is an interesting new phenomenon that, like triangular 
flow, results from initial state fluctuations and subse- 
quent hydrodynamic evolution, and provides the last puz- 
zle piece needed for a complete description of long-range 
dihadron correlations. Unlike the usual directed flow, 
this phenomenon has no correlation with the reaction 
plane and should depend weakly on rapidity. It changes 
sign around p t ~ 0.8 GeV, and may reach values as large 
as 10% at high p t for mid-central collisions. 

We have proposed a specific method to analyze this 
new observable, which eliminates correlations from mo- 
mentum conservation. Such dedicated measurements of 
this new v\ at RHIC and LHC will help constrain mod- 
els of initial conditions and confirm the hydrodynamic 
behavior of the collision system. 
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1 Note that i> EP>1 , does not actually denote a plane, but is a emphasize the similarity to a v 2 event plane analysis 

specific direction for directed flow. The terminology is meant to 



